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Abstract 

We study the speed of convergence of the explicit and implicit space-time discretization 
schemes of the solution u(t, x) to a parabolic partial differential equation in any dimension 
perturbed by a space-correlated Gaussian noise. The coefficients only depend on u(t, x) 
and the influence of the correlation on the speed is observed. 
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1 Introduction 

Discretization schemes for parabolic SPDEs driven by the space-time white noise have been 
considered by several authors. I. Gyongy and D. Nualart [2] and [10], have studied implicit time 
discretization schemes for the heat equation in dimension 1. J. Printems [15j has studied several 
time discretization schemes (implicit and explicit Euler schemes as well as the Crank-Nicholson 
one) for Hilbert- valued parabolic SPDEs, such as the Burgers equation on [0,1], introduced 
several notions of order of convergence in order to deal with coefficients with polynomial growth 
and proved convergence in the Hilbert space norm. This work has been completed by E. 
Hausenblas [TT] . who studied several schemes for quasi-linear equations driven by a nuclear 
noise, and taking values in a Hilbert or a Banach space X. Several approximation procedures 
(such as the Galerkin approximation, finite difference methods or wavelets approximations) were 
considered, but the coefficients of the SPDE were supposed to depend on the whole function 
u(t, .) in X, and not only on (t,x). Notice that, unlike [TT], the coefficients considered in this 
paper do not depend on the whole function u(s, .). 

I. Gyongy [7j has studied the strong speed of convergence in the norm of uniform convergence 
over the space variable for a space finite-difference scheme u n with mesh 1/n for the parabolic 
SPDE with homogeneous Diri chief's boundary conditions. He has also studied the speed of 
convergence of an implicit (resp. explicit) finite-difference discretization scheme u n,m (resp. 



u^) with time mesh T/m and space mesh 1/n for the solution u to the following parabolic 
SPDE in dimension 1 driven by the space-time white noise W: 

§ (t, x) = 0(t, x) + a(t, x, u(t, x))gg + b(t, x, u(t, x)) , , s 

u[t,0) =u(t,l) = 0, 1 ' J 

with the initial condition uq. He has proved that, if the coefficients a(t, x, .) and b(t, x, .) satisfy 
the usual Lipschitz property uniformly in (t, x) and if the functions a(t, x, y) and b(t, x, y) are 
1/4-Holder continuous in t and l/2-H61der continuous in x uniformly with respect to the other 
variables, then for t G]0, T], p G [1, +oo[, < (3 < \ and < 7 < ~ one has: 

sup E(|w n ' m (t,x) -M(t,x)| p ) <K(t) (m-^ + n-^) . (1.2) 

a:G[0,l] 

Furthermore, if u G C 3 ([0, 1]), then (II. 2p holds on [0, T] with /3 = |, 7 = | and with a constant 
X which does not depend on t. A similar result holds for the explicit scheme u 1 ^ if < q < |. 

A. Debussche and J. Printems |5j have implemented simulations of a discretization scheme 
for the KDV equation, and C. Cardon- Weber [2] has studied explicit and implicit discretization 
schemes for the function-valued solution to the stochastic Cahn-Hilliard equation in dimension 
d < 3 when the driving noise is the space-time white noise. The polynomial growth of the drift 
term made her require the diffusion coefficient a to be bounded, and she proved convergence 
in probability (respectively in L p with a given rate of a localized version) of the scheme. 

In the present paper, we deal with a rf-dimensional version of (11.11) . As it is well-known, we 
can no longer use the space-time white noise for the perturbation; indeed, in dimension d > 2, 
the Green function associated with — A with the homogeneous Dirichlet boundary conditions 
011 [0, l] d is not square integrable. Thus, we replace W by some Gaussian process F which is 
white in time and has a space correlation given by a Riesz potential f(r) = r~ a , i.e., such that if 
A&ndB&re bounded Borel subsets of R d , E(F(s, A) F(t, B)) = (s At) J A dx f B dy\x-y\' a for 
some a G]0,2 A d[. See e.g. [T2], [1], [TJj and [3] for more general results concerning necessary 
and sufficient conditions on the covariance of the Gaussian noise F ensuring the existence of a 
function-valued solution to (11. If) with F instead of W . 

The aim of this paper is threefold. We at first study the speed of convergence of space 
and space-time finite discretization implicit (resp. explicit) schemes in dimension d > 1, i.e., 
on the grid (g, , 1 < k < d)), < i < m, < j k < n and extended to [0,T] x [0, l} d 
by linear interpolation. As in [7] and [8], the processes u n and u n,m (resp. w^J have an 
evolution formulation written in terms of approximations {Gd) n , {Gd) n,m and {Gd)m of the 
Green function Gd, while u is solution of an evolution equation defined in terms of Gd- These 
evolution equations involve stochastic integrals with respect to the worthy martingale-measure 
defined by F (see e.g. [18J and [4|). As usual, the speed of convergence is given by the norm 
of the differences of stochastic integrals; more precisely, the optimal speed of convergence for 
the implicit scheme is the norm of the difference Ga(.,x,.) — (Gd) n ' m (-, x, .) in L 2 ([0, T], Hd), 
where Tid is the Reproducing Kernel Hilbert Space defined by the covariance function. More 
precisely, if and ip are continuous functions on Q = [0, l] d , set 



<¥^>n d = \ \ p(x)f(\x-y\) ip(y)dxdy. (1.3) 

J Q JQ 

We denote by TCd the completion of this pre-Hilbert space; note that Tid elements which are not 
functions and that a function ip belongs to TCd if and only i£ Jg Jg\(p(y)\ f (\y — z\) \<p(z) \ dy dz < 
+00. However, unlike in [7j and [8], the functions <fj(x) = V2 sin(j7rx), j > 1 and (fj(K n (x)), 
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1 < j < n, where K n (y) = [ny] rT l are not an orthonormal family of 7i\. Thus, even in 
dimension d = 1, the use of the Parseval identity has to be replaced by more technical com- 
putations based on Abel's summation method. Similar results could be obtained for a more 
general space covariance, provided that it is absolutely continuous and that its density / sat- 
isfies some integrability property at the origin (see e.g. [I], [E]). However, the speed of 
convergence would depend on integrals including /, which would make the results less trans- 
parent than that stated in the case of Riesz potentials. The key technical lemmas, giving upper 
estimates of \\G d (.,x, .) - (G d ) n {. } x, .)||l9([o,do[,?^) and || (G d ) n (., x, .) - (G d ) n ' m (., x, .)\\m[o,T],H d ) 
(resp. \\(G d ) n (.,x, .) - (G d ) n m (. } x, .) \\^([o,T],H d )), are proved in section H 

We describe the discretization schemes in any dimension d > 1 and introduce some notations 
in section 2. In section [31 an argument similar to that in [7] shows that for < a < d A 2, and 
p 6 [1, +oo[, if m is regular enough, then 

sup E(\\u(t,x) -u n (t,x)\\ 2p ) < C p , a n- {2 ~ a)p , (1.4) 

(t,x)£[0,+oo[xQ 

and extending [S] we prove in section 4 that 

sup E{\\u n {t,x) -u n ' m {t,x)\\ 2p ) < C p>a m-( 1 -^ p . (1.5) 

(t,x)e[0,T]xQ 

If d = 1, as a y 1 the space density becomes more and more degenerate and the speed of 
convergence approaches that obtained by Gyongy for the space-time white noise. 

In dimension d > 2, the proof depends on the product form of the Green function and its 
approximations, as well as of upper estimates of \x — y\~ a in terms of Yli=i \ x i ~ Vi\~ ai f° r some 
well-chosen aij. Thus, estimates of the Hd-norm of the differences of G d (s,x, .) — (G d ) n (s, x, .), 
(G d ) n (s,x, .) — (G d ) n ' m (s, x, .) and (G d ) n (s,x, .) — (G rf )^(s,x, .) in dimension d > 2 depend on 
bounds of the Hi-norm of similar differences as well as of 7i r -norms of G(s, x, .), G n (s, x, .) and 
G n ' m (s, x, .) for r < d. 

Section [5] contains some numerical results. For T = 1, we have implemented in C the (more 
stable) implicit discretization scheme for affine coefficients a(t, y,u) = a\u + a 2 and b(t, x, u) = 
bi u + 62 and for a(t,y,u) = b(t,y,u) = a + b cos(w). We have studied the "experimental" 
speed of convergence with respect to one mesh, when the other one is fixed and gives a "much 
smaller" theoretical error. The second moments are computed by Monte-Carlo approximations. 
These implementations have been done in dimension d — 1 for the space-time white noise W 
and the colored noise F. As expected, the observed speeds are better than the theoretical ones, 
and decrease with a. For example, choosing and M "large" with M > N 2 and considering 
"small" divisors n of N, we have computed the observed linear regression coefficient and drawn 
the curves of sup xg [ 01 ] hx(E(\u n,M (1, x) — u N,M (l,x)\ 2 )) as a function of ln(n) for various values 
of a. 

Note that all the results of this paper remain true if in (11. lj) we replace the homogeneous 
Dirichlet boundary conditions u(t,x) = for x G dQ by the homogeneous Neumann ones 
|^(t,x) = for x G dQ. In this last case, the eigenfunctions of — A in dimension one is 
(fo(x) = 1 and for j > 1, (fj(x) = \pl cos(j7rx). Since the upper estimates of the partial sums 
^2f = i>-Pj(x) used in the Abel transforms still hold in the case of Neumann's conditions, the 
crucial result is proved in a similar way in this case, and the speed of convergence is preserved. 
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2 Formulation of the problem 

Let (Q, J 7 , P) be a probability space, Q = [0, l] d for some integer d > 1 and let F = (F((p) , 9? G 
P(R + x Q)) be an L 2 (P)-valued centered Gaussian process, which is white in time but has a 
space correlation defined as follows: given ip and ip in D(R + x Q), the covariance functional of 
F{ip) and F(i/j) is 

J(cp,1>)=E(F(<p)F(4>))= dt / ip(t,y)f(y-z)^(t,z)dydz, (2.1) 

JO J J{Q-Q)* 

where (Q — Q)* = {y — z : y,z G Q, y 7^ z} and f : (Q — Q)* — > [0, +oo[ is a continuous 
function. The bilinear form J defined by (12.1 ft is non-negative definite if and only if / is the 
Fourier transform of a non-negative tempered distribution \x on Q. Then F defines a martingale- 
measure (still denoted by F), which allows to use stochastic integrals (see [IE])- In the sequel, 
we suppose that for ^ 0, f(z) = \z\ a , where \z\ denotes the Euclidean norm of the 

vector z. Since x 2 + y 2 > 2xy, if 01 j = a2 _jf for 1 < j < d and ad = a2~ d+1 , there exists a 
positive constant C such that for any z — (zx, ■ • ■ , Zd) G R d , 

f{z)<C J] f a .( Zj ), (2.2) 

l<3<d 

where / a (C) — |C| _Q f° r an y C G ^ C 7^ 0- To lighten the notations, for this choice of / and 
ip EHd set 

IMI(«)=/ / \<p(y)\\y-z\- a \<p(z)\dydz. (2.3) 

For any t > 0, we denote by jF t the sigma-algebra generated by {P([0, s] x A) : < s < t , A G 
Q}. Let a : [0, +oo[xQ x R — > R and 6 : [0, +oo[xQ x R — > R. Suppose that there exists a 
positive constant C such that for s,t G [0, 00 [, x, ?/ G Q, r, t> G R, the linear growth condition 



(E3D and either Lipschitz condition 023]) . (1231) or (EZD hold 

|a(t,x,r)| + |6(t,x,r)| < C (1 + |r|) , (2.4) 

and for D(s, t, x, y, r, v) = \cr(s, x, r) — <r(t, y,v) \ + \b(s, x, r) — b(t, y, v) \ 

D(t,t,x,x,r,v) <C\r-v\, (2.5) 

D(t,t,x,y,r,v) < C(|x - y| 1_ § + \r-v\), (2.6) 

D(s,t,x,y,r,v) < C(\t — s\^~^ + \x — y] 1 ^^ + \r — v\). (2.7) 



For any function u$ which vanishes on the boundary of Q, let u(t, x) denote the solution to the 
parabolic SPDE, which is similar to (11. ip 

f §jf (*, x) = Au(t, x) + cr(t, x, u(t, x))j££ + b{t, x, u(t, x)) , , . 

\ u(t, x) = for x G dQ , 1 ' ' 

with initial condition u(0,x) = uq(x). Let N* denote the set of strictly positive integers. For 
any j G N* and £ G R, set <^-(£) = V"2sin(j7r£) and for = (/ci, • • • , G N* d , set 

d d 

\k\ = ^2kj, fk (x) = Yl Vkjfa.) for x = (xi, . . . ,x d ) G R d . 
i=i i=i 
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Let Gd{t, x, y) denote the Green function associated with J^— A on Q and homogeneous Dirichlet 



boundary conditions; then for t > 0, x,y G Q, Gd(t,x,y) 
and 



\G d (t,x,y)\ <Ct 2 exp ( - c 



\x 



J2keN* d ex P( 

2 



y 



t 



-)• 



|fc| 2 7r 2 t)^fc(a;)^fc(|/) 



(2.9) 



When d = 1, set G\ = G. These upper estimates are classical when the domain Q has a smooth 
boundary under either homogeneous Neumann's or Dirichlet's boundary conditions (see e.g. [6], 
[13J). A simple argument shows that they can be extended for these homogeneous conditions 
on the set Q = [0, l] d ; see e.g. pQ for the similar case of the parabolic operator + A 2 on 
[0, 7r] d . The equation (12.81) makes sense in the following evolution formulation (see e.g. [IB] for 
(2=1): 



u(t,x)= / G d (t,x,y)u (y)dy 



G d (t - s,x,y) 



x [<r(s, y, u(s, y)) F(ds, dy) + b(s, y, u(s, y)) dsdy] 



(2.10) 



We also consider the parabolic SPDE with the homogeneous boundary conditions |^(t, x) = 
for x G dQ. Then the functions (cpj ; j > 1) are replaced by </?o(0 = 1 and (fj(0 = \/2 cos(j7r£) 
for (gI and j > 1. All the other formulations remain true with k G N rf instead on W d . 



2.1 Space discretization scheme 

As in [TJ, we at first consider a finite space discretization scheme, replacing the Laplacian by 
its discretization on the grid = = (— , • • • , where kj G {0, • • • , n}, 1 < j < d. In dimension 
1, we proceed as in [7J, and consider the (n — 1) x (n — l)-matrix D n associated with the 
homogeneous Dirichlet boundary conditions and defined by D n {i,i) = —2, D n (i,j) = 1 if 



1 and D r 



for 



j\ > 2; then 9 is replaced by n 2 D n u n (t, .), where 



u a (t) denotes the (n — l)-dimensional vector of an approximate solution defined on the grid 
j/n, I < j <n. In arbitrary dimension, we proceed as in [2] and define Dn by induction. Let 

ma trix. Let 
matrix A, let diag(A) 



D n l) = D n and suppose that D n d has been defined as a (n 



\d-i 



x n 



1) 



x 

2-1 



7(4 denotes the k x k identity matrix and given a (n — 1 

denote the (n — l) d x (n — l) d matrix with d — 1 diagonal blocs equal to A; let D„ denote the 



n 



l) d x (n — l) d -matrix D« defined by 



/ — 1Id n d-\ Id r 



D«> = diag{D n ^ 



(d-ih 



V 



Id n d-1 







-1Id n d-\ 



Id n d-1 



Id n d-1 







~1Id n d- 
Id„d-1 



Id n d-\ 
-II d n d-l 



J 



Mn(i,Xk). witn *k 



n 
[n- 



Let vP'it) denote the (n — l) d dimensional vector defined by vP'it)^ 

(xk-L, ■ ■ ■ ,Xk d ), where kj is the unique integer such that x^ = and kj G {1, • • • 
such that k = (fc d -l)(n-l) d_1 + - • ■ + (k 2 -l){n-l) + k 1 . Let £ = {x k : k G {1, 
□x k be the lattice parallepiped of diagonal x k = (x^, ■ ■ • , Xk d ) and (xfc 1 + -,••• , Xk d + -), and 
set~F n (t,Xk) = Ja x dF(t,x). Given a function /i : [0,+oo[xQ x R -> R, and w G M r , let 
/i(t, x, tt) = (/i(t, x, Mi), • ■ ■ , h(t, x, u r )). Then vP'it) is solution to the following equation 



- 1} is 

in}, 



diT{t) = n 2 D n d) tT{t) dt + na(t, x, tT{t)) dF(t, .), +b(t, x, iT(t)), 



(2.11) 
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u" 



(0) = (uq (-) , j G {1, • • • , n — l} d ) . We then complete u n (t, .) from the lattice C to Q as 
follows. If d = 1, set w n (t,0) = u n (t,l) = 0, w n (t, ^) = k^O) = [ny)/n, <p?(i/n) = 

<Pj(i/n) for < i < n, and for x &]i/n, (i + l)/n[, < i < n, 1 < j < n — 1, let y?™(x) = 
^ (i) + (nx - i) h - <Pj (£)] , and let 



A' 



-4 sin — 

2n 



n 



j 2 7t 2 4 with 4 



denote the eigenvalues of n 2 D n = n 2 Dn^', then for t > 0, x, y G [0, 1] 



■ 2 F 
S1U ' 2^ 



J7T 

2n 



7T 



n-1 



(Gi) B (t, x, y) = £ exp(A™ t)^(x) ^K(y)) . 



(2.12) 



In dimension d > 2, we also complete the solution u n (t,x) from x G £, defined as u n (t,x k ) = 
to x G Q by linear interpolation, interpolating inductively on the points (x, y) for x G R 
and y = (k i+ i/n, ■ ■ ■ , kd/n). The eigenvalues and eigenvectors of n 2 D„ are A£ = X)j=i ^% 
and <^ k ftf , ■ ■ • , ) • For t > 0, x and y G Q if K n {y) = (/^(j/i), • ■ ■ , k„(^)), let 



{G d ) n {t,x,y) = exp(A^K(xK(K„(y))- 

ke{l,»- ,(n-l) d } 



(2.13) 



When d = 1, simply set Gy — G and (Gi) n = G n . Then the linear interpolation of u n (t, .) from 
the lattice £ to Q = [0, l] d is solution to the evolution equation 



u n (t,x) = / (G d ) n (t,x,y)u {K n (y))dy + 



Q 



(G d ) n (t-s,x,y) 



x 



a(s, n n (y),u n (s, K n (y))F(ds, dy) + b(s, K n (y),u n (s, K n (y)) dsdy 



(2.14) 



The n x n matrix D n = D n associated with the homogeneous Neumann boundary conditions 



-1, D n (l,2) 



D n (n,n - 1) 



1 and for 2 < i < n — 1 
i\ > 2. 



is defined by D n (l, 1) = D n (n,n) 

and 1 < j < n, D n (i,i) = -2, D n (i,j) = 1 if |j - i| = 1 and D n (i,j) = for |j 
The inductive procedure used to construct Dn is similar to the previous one, replacing 1 by 
Id n a. Then the eigenvalues of n 2 D n are A" = — 4n 2 sin 2 (||) = — j 2 it 2 c£ with c J n G [^s, 1]- 
The corresponding normed eigenvectors (ej, < j < n — 1) are again evaluations of tpj. More 
precisely, ej(fc) = 77^ ^'(^f^O for < j < n — 1 and 1 < < n. The eigenvalues A- and the 

eigenf unctions </2k of n 2 are defined in a way similar to the Dirichlet case, taking sums over 
kG {!,■■■ , n d }. Formulas similar to (12.12ft and (12.131) still hold and (12. 141) is unchanged. 



2.2 Implicit space-time discretization scheme 

We now introduce a space-time discretization scheme. Given T > 0, n, m > 1 we use the space 
mesh 1/n and the time mesh T/m, set tj = iTm~ l for < i < m and replace the time derivative 
by a backward difference. Thus for d = 1, in the case of Dirichlet's homogeneous boundary 
conditions, set wo = (v>o{j/ n )> 1 <J <n — 1) and for i < m, set ttj = {u n,m {iTm~ l , jn _1 ) , 1 < 
j < n- 1), and for g = a and p = b let g-(^, .,«») = g^UJu" 1 , {u n ' m {t h jn' 1 )) , 1 < j < n-1). 
Let \3 n m F(ti, ■) denote the (n — l)-dimensional Gaussian vector of space-time increments of F 
on the space-time grid, i.e., for 1 < j < n — 1, set 

n nim F{t h j) = nmT-^FiU+u (j + l)^ 1 ) - F{t h {j + l)^ 1 ) - F(t t+1 , jn" 1 ) + jn' 1 )}; 
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then for every < i < m 



T 

u i+ i = Ui + n 2 — D n u i+1 + TmT 1 [a{% .,Ui)n n , m F(t h .) + b(U, .,«»)] . (2.15) 



Since Id — Tm D„ is invertible 



rj~t 1 rj~t 

u i+1 =(Id - -D n )-^u + J2(Id - -D n )-«- k -V [a(t k , ., u k p n>m F(t k , .) + b(t k , ., u k )} 

k=0 



(2.16) 



If d > 2, we set □„, m -F(t i , x k ) = n d mT 1 f 1+1 J n dF(t,x), and for homogeneous Dirichlet's 

(resp. Neumann's) boundary conditions, define similarly u i+ i as the (n — l) d -dimensional (resp. 
n d -dimensional) vector such that (I2.16P holds with instead of D n . We only describe the 
scheme in the case of Dirichlet's conditions; the case of Neumann's conditions is obviously dealt 
with by obvious changes. The process u n ' m is defined on the space-time lattice Ct = {(^,x k ) : 
< i < m, k G {1,- •■ ,(n - l) d }} as (u n ' m (t h x k ) , < i < m , k G {1, • • • , (n - l) d }) = ^ ; 
it is then extended to the time lattice (U, x), < i < m, x G Q as in the previous subsection, 
and then extended to [0, T] x Q by time linear interpolation. Since Ak = J2i=i ^k anc ^ ^(xjj 



are the eigenvalues and eigenvectors of , if 



(G d ) n ' m (t } x,y) = i 1 - Tm~ l K)- [mtT ~ 1 <pl(x) ^ n {y)) , (2.17) 

ke{l,- ,(n-l) d } 

then for t = iTm' 1 , 1 < i < m, if for s G [0,T], one sets A m (s) = [msT^ 1 ]m^ 1 one has: 

u n > m (t,x)= [ (G d r m (t,x,y)uo(K n (y))dy+ f [ (G d ) n,rn (t - s + —,x,y) 
Jq Jo Jq m 

x cr(A (y)))F(ds,dy) + b(A (y)))dyds\. 

(2.18) 

Again for d = 1, let G n ' m = (d 



\n,m 



2.3 Explicit schemes 

For T > 0, a space mesh n~ x and a time mesh Tm -1 , we now replace the time derivative by a 
forward difference. Thus if u 1 ^ denotes the approximating process defined for t = U = iTmT 1 
and x kj G {1, • • • , n — 1}, setting = u^fa, .), we have 

u l+l = ^ + n 2 Tm~ 1 D { ^u i + TrrT 1 [a(t h ., r^n^Ffa, .) + b(t h .,u { )} . (2.19) 

In the case of homogeneous Dirichlet boundary conditions, let (G^^it, x, y) denote the corre- 
sponding approximation of the Green function G d defined by 

(G d r m (t,x,y) = (1 + Tm- 1 A k ) MT_1] <p%x) ^ n {v))- (2.20) 

ke{i,- ,(n-l) d } 

Again for d — 1, let G 1 ^ = (Gi)^. Then for t = ti = iTm~ l , when completing the solution 
w^(tj, .) from the space lattice £ to Q, we obtain the solution to the following equation 

<(t,x)= / (G d )l(t,x,y)u (K n (y))dy+ [ [ {G d )l(t - s + T/m,x,y) 



Q JO 
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cr(A (y)))F(ds,dy) + b(A (y)))dyds 

(2.21) 

We then complete the process by time linear interpolation and obvious changes yield 

the explicit scheme for homogeneous Neumann boundary conditions. 



3 Convergence results for the discretization schemes 

In this section, we study the speed of convergence for the <i-dimensional space scheme and then 
of the <i-dimensional implicit and explicit space-time schemes. For the sake of simplicity, we 
only write the proofs in the case of homogeneous Dirichlet boundary conditions. The following 
result states that the solutions u, u n , u n,m and u 7 ^ exist and have bounded moments uniformly 
in n,m. The proofs for u can be found in [T3]; see also [I] and [3]. The arguments for the 
approximations are similar using (1A.9[) . (lA.lOp . flA.13[) and OA. 140 and the version of Gronwall's 
lemma in stated in [H] Lemma 3.4. 

Proposition 3.1 Let u$ G C(Q) satisfy the homogeneous Neumann or Dirichlet boundary con- 
ditions, and suppose that the coefficients a and b satisfy the conditions $2.4]) and $2. 5]) ; then 
the equation $2.10]) (resp. $2. 14 ), $2.181) and $2.21]) ) has a unique solution u (resp. u n , u n ' m 
and vJ^J such that for every pG [1, +oo[ and T > 0: 

sup sup sup E(\u{t,x)\ 2p + \u n (t,x)\ 2p + \u n > m (t,x)\ 2p ) + \u n m {t,x)\ 2p ) < oo . (3.1) 

n,m 0<t<T xeQ 

We now prove Holder regularity properties of the trajectories of u and u n . Note that for 
u, a similar result has been proved in [T7j for the heat equation with free boundary with a 
perturbation driven by a Gaussian process with a more general space covariance; see also [3] 
for a related result in the case of a more general even order differential operator. 

Proposition 3.2 Suppose that the coefficients b and o satisfy the Lip schitz property $2. 5\) . that 
the initial condition Uq satisfies the homogeneous Dirichlet or Neumann boundary condition. 

(i) Suppose furthermore that Uq G C l ~?(Q) and fix T > 0. Then, for every p e [l,+oo[, 
there exists a constant C such that for x, x' G Q and < t < t' <T, 

sup E (\u(t, x) - u(t, x')\ 2p ) <C\x' - x\ p{2 ~ a) , (3.2) 

0<t<T 

supE(\u{t',x) -u(t,x)\ 2p ) <C\t' -t\ p{1 ~^ . (3.3) 

x&Q 

(ii) Suppose furthermore thatu G C 2 (Q); then for every p G [1, +oo[, there exists a constant 
C such that for x, x' G Q and < t < t! < T, 

supsupE(|n n (t',a;) -u n (t,x)\ 2p ) <C\t'-t\ p ^ (3.4) 

n>l xSQ 

sup sup E(\u n (t,x') -u n (t,x)\ 2p ) < C\x' - x\ p{2 ~ a) (3.5) 

n>l 0<t<T 

Proof: The proofs of (13.21) and (13.31) can be adapted from Sanz-Sarra [T7j (see also [S]), and 
are therefore omitted. For the sake of completeness, we sketch the proof of (13.41) . For every 



S 



t > 0, let v n (t, x) = J (Gd) n (t, x, y) u (K n (y) dy and w n (t, x) = u n (t, x) — v n (t, x), where (G d ) n 



is the fundamental solution of ^ — A n = 0, 



A„u ( y) ^ £ [u( £ Me 3 + M±ie.) _ ^(M) + „(£ M c , + Mzi 

(3.6) 



n n 



and (e^, 1 < i < d) denotes the canonical basis of ¥L d . Then if uq G C 2 (Q), 

v n (t,x)=u (x)+ [ [ (G d ) n (s,x,y)A n u (y)dyds. (3.7) 

Using the fact that A n u is bounded if u G C 2 (Q), and (IA.22[) . we deduce that for any A > 

sup sup \v n (t, x) — v n (t', x)\ < C \t' — . (3.8) 

Computations similar to those used in [17], using Burkholder's and Holder's inequalities with 
respect to suitable measures, flA~22l - (lA~23l) and (EH]) , show the existence of C p > such that 
for any < t < t' < T 

supsupE(K(t,x) -w n {t',x)\ 2p ) < C p \t' - t\ p{1 ~^ . (3.9) 

n>l xgQ 

The inequalities (13. 8p and (13.91) conclude the proof of (13 .4p . □ 
The first convergence result of this section is that of u n to u. 

Theorem 3.3 Let o and b satisfy the conditions \2.1$ and ^2. 61) . u and u n be the solutions to 
Ii2.10\) and \2.1J$ respectively, where the Green functions G d and {G d ) n are defined with the 
homogeneous Neumann or Dirichlet boundary conditions on Q. 

(i) If the initial condition Uq belongs to C 3 (Q), then for every T > and p G [1 + oof, there 
exists a constant C P (T) > such that: 

sup E(\u{t,x) -u n {t,x)\ 2p \) < C p {T)n^ 2 - a)p . (3.10) 

(t,x)e[0,T]xQ 

(ii) If the initial condition uq belongs to C 1_ ^(<5), then there exists v > such that given any 
p G [1, +oo[, there exists a constant C p > such that, for every t > 0: 

sup E(\u(t,x) -u n {t,x)\ 2p \) <C p r v n -^ p . (3.11) 

x<=Q 

(Hi) Finally, if uq belongs to Cq(Q), then for all p G [1, +oo[ ; as n — > +oo, sup E,[\u(t,x) 

{t,x)e[0,T]xQ 

— u n {t, x)\ 2p \) converges to 0, and the sequence u n {t,x) converges a.s. to u(t,x) uniformly on 
[0,T]xQ. 

Proof: As in [7], set v(t,x) = J Q G d (t,x,y)u (y) dy, v n (t,x) = j Q (G d ) n (t,x,y) u (n n (y)) dy, 
and u{t,x) = v{t,x) + w(t,x) u n {t,x) = v n (t,x) +w n (t,x). If u G C l ~?{Q) (and hence is 
bounded), using ( 14. ip and ( 1A.9[) . we deduce that for any A g]0, 1[, there exists /i > 0, C > 
such that for t > 0, v = A V //, 

supK*,x) -v n {t,x)\ < / \\G d {t,x,y) - {G d ) n {t,x,y)\ \u Q {y)\ 

x£Q JQ l 
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+ \(G d )%x,y)\\u (y) -u (K n (y))\ dy < C (1 + t~ v ) e^n^f ) 



(3.12) 



If uo e C 3 (Q), then since (resp. (G d ) n ) is the fundamental solution of — A = (resp 
— A n = 0), where A n is defined by ( 13.61) . integrating by parts we deduce that v(t,x) = 
u (x) + $1 f Q G d (s, x, y) Au (y) dy and v n (t, x) = u (x) + f Q J Q (G d ) n (s, x, y) A n u (y) dy. Hence 

\v(t,x) — V n (t,x)\ < Y^=l-Ai(t' X )i wnere 



Ax(t,x) = \uo(t,x) - uo(«n(a;))| , A 2 (t,x) 




[G d (s,x,y) - (G d ) n (s,x,y)}Au (y)dyds 



JQ 



A 3 (t,x) 




(G d ) n (s,x,y)[Au (y) - A n u Q (K n (y)))dyds 



o JQ 



< Cn- 1 , the 



Since Au Q i s bo unded and \\u (.) - m (k„(.)) ||oo + ||Aw (.) - A n u (/c n (.)) 
inequalities (1A.9|) and (14. 2 p imply 

sup \v(t, x) — v n (t, x)\ < C rT l . 

(t,x)e[0,+oo[xQ 

Furthermore, for every t E]0,T], sup^Q E(|u>(t, x) — w n (t,aOP*)<C£° =1 fli(t), where 



(3.13) 



Bx(t) = supE( 
B 2 (t) = supE 



B 3 (t) =supE 
5 4 (t) =supE 
B 5 (t) =supE 

x&Q 

B a (t) =supE 




G d (t - s,x,y)(a(s,y,u(s,y)) - a(s, n n (y),u(s, K n {y)))F(ds,dy) 



2p 



OJQ 



G d (t - s, x, y)(cr(s, K n (y), u(s, n„{y))) 



\ 2p 

- a(s, K n (y), u n (s, K n (y))JF(ds, dy) 

G d (t - s, x, y) - {G d ) n (t - s, x, y))a(s, K n (y), u n (s, K n (y))F(ds, dy) 
G d (t - s, x, y) b(s, y, u(s, y)) - b(s, K n (y), u(s, K n (y)) )dyds 



2p 




G d {t - s,x,y)[b(s, K n (y),u(s, K n (y))) - b(s, n n (y),u n (s, K n (y)) )dyds 

2p 



2p 



G d (t - s, x, y) - (G d ) n (t - s, x, y) )b(s, n n (y),u n (s 1 n n (y))dyds 



(3.14) 



Burkholder's inequality, (IA. 11) . Holder's inequality with respect to \\G d (t — s, x, .)||? a ) ds, Fubini's 
theorem, ( 12.61) . Schwarz's inequalities and ( 13. 2ft imply that 



B^t) < C p supE / \G d (t-s,x,.)\(n-— + \u(s,.) -u(s,K n (.)\) 

x&Q 



ds 



(a) 



<C P (t- sY 



n 



-p(2-a) 



+ SUp E(|li(s,f) - u(s,Kn(g))\ 



2p\ 



(s,f)£[0,i]xQ 



ds < C n n^ 2 - a) . 



(3.15) 



Similar arguments based on (lA.lft and (12.61) (resp. (14.31) . (12.41) and (13. ip ) imply that 

B 2 (t) <C P [ (t-s)-f (swp\v(s,x) -v n (s,x)\ 2p + supE(\w(s,x) - w n ( Sl x)\ 2p ))ds, (3.16) 
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B 3 (t) <C p n-V-°» p . 



(3.17) 



The deterministic integrals are easier to deal with. Using Holder's inequality with respect to 
the measure \G(t — s,x,y)\ dyds, (12.91) . (12.61) and (13. 2ft we deduce that 



B 4 (t) <C P sup e( I i \G d (t - s,x,y)\dyds) [ [ \G d (t - s,x,y)\ 
x<eq v Jo Jq ' Jo Jq 

x (n~ {2 ~ a)p + E(\u(s,y) -u(s, K n {y)) \ 2p )}dydzds < Cn" p{ - 2 ~ a) . 
Similarly, Hdlder's inequality, (JUS]) and flUD (resp. (J23D, flUD and g^D) imply 

ds, 







f 


sup 


lo 


- x&Q 



\2p 



\2p 



x&Q 



B 6 (t) <Cn~ 2p . 



(3.18) 

(3.19) 
(3.20) 



The inequalities (13. 151) - (13. 201) imply that for any T > and p G [1, +oo[, there exists a constant 
C > such that for < t < T, 



supE(Kt,x) -w n (t,x)\ 2p ) < C[n- p{2 - a) + / (t-s)~% 

x&Q ^ 







X 



sup \v(s,x) - v n (s, x)\ 2p + supE(|w(s, x) - w n (s,x)\ 2p ) 

x£Q x£Q 



ds). 



(3.21) 



Thus, (I3.13P and Gronwall's lemma (see e.g. [8], lemma 3.4) imply that if uq G C 3 (Q), 



supE{\w(t,x) -w n (t,x)\ 2p ) < C p 



n 



-p(2-a) 



+ / (t-s) 



x sup E( | it; (s,x) - u> n (s, x)\ 2p ) ds 



a;<E( 



< C p n- p{2 - a) . 



This inequality together with (l3~T3l yield (l3~T0l) . If w G C 1- i (<2), using again Gronwall's lemma 
and (13.121) . we deduce that for some A g]0, 1[, one has 



supE{\w{t,x) -w n {t,x)\ 2p ) <Cp(n- pi2 - a ~> + \ (t - s)"t 







s -A n -p(2-a) 



+ supE(|u>(s,x) -w n (s,x)| 2p ) 



This inequality and (13. 1 2[) imply (13. lip . 

Finally, let u G C°(Q) and for any e > 0, let u 0i£ denote a function in C 3 (Q) such that 
11% — M o,e||oo < £• Let m £ = f £ + w e and = t> ™ + w™ denote the previous decompositions of 
the solution u £ and its space discretization with the initial condition uo, £ - Then 



sup \v (t, x) — v n (t, x)\ < sup \v £ (t, x) — v™(t,x)\ 

(t,x)e[0,T]xQ (t,x)e[0,T]xQ 



+ 



G d (t, x, y) \u (y) - u Jy))dy 



<Ce + sup \v e (t, x) — v™(t, x)\. 

{t,x)e[0,T]xQ 



(G d ) n (t, x, y) [u (K n (y)) - u 0)£ (K n (y))]dy 

(3.22) 
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Hence (15211) and 113221) imply that 



supE(|w(t,x) -w n (t,x)| 2p ) < C 



£ + n' 



-p(2-a) 



+ / (t-s)-^swpE(\w(s,x) -w n (t,x)\ 2p )ds 



Gronwall's lemma concludes the proof of the theorem. □ 

We now prove the convergence of u n ' m and of nj^ to u n as m — > +00. 

Theorem 3.4 Lei cr and b satisfy the conditions pO j and (IT7). Taen 

// Mo ^ C 2 (<5) ; ^ en / or ever U T > and p G [1, +00 [, iaere exists a constant C P (T) > 
such that 

sup sup supE(\u n (t,x)-u n > m (t,x)\ 2p ) <C p (T)m- p{1 -^ . (3.23) 
i>i te[o,T] xeQ 

(ii) Ifuo G C(Q), then sup n>1 sup tg r 0T i sup^g |n n (t,x)— n n,m (t,x)| converges to as m — > +00 
and /or even/ 1 > and p G [1, +oo[ there exists a constant C p (t) such that 

sup supE(|n n (t,x) -n n ' m (t,x)| 2p ) < C p (t) m- p(1 ~f > . 

n>l xSQ 

(m^ T/ie results of (i) and (ii) hold with u™ t instead of u n ' m if one requires that < q < ~. 

Proof: Let v"(f,a;) = f Q (G d ) n (t, x, y)u (K n (y))dy, v n > m (t,x) = f Q (G d p m (t, x,y)u (K n (y))dy. 
Suppose at first that no G C 2 (<3) and as in the proof of (3.23) in [8], for d = 1 set / = 
sup te[0 ,T] sup^y |f n ' m (t,x) - f n (t,x)| < J2 3 i= i h, where 





= sup 


sup 






te[o,T] 


se[o,i] 




h 


= sup 


sup 






te[o,T] 


X6[0,l] 




h 


= sup 


sup 






te[o,T] 


ze[o,i] 





The inequalities (3.27) and (3.28) in jH] imply that I2 + I3 < C m 2. Furthermore, using an 
estimate of [8], we deduce that 



n-1 



J x <C sup sup sup J~]j 2 exp (A™ [^1 

n>l te[0,T] x6[0,l] ~i V 



T m 



1 — exp 



I 1 \ J m V J m 



For t <Tm 1 , [^] = and the right hand-side of the previous inequality is 0. If t > Tm 1 . 
then there exists a constant c > such that ^ [^r] > ct and using (1A.7j) we deduce that 



n-1 



Ji <C sup sup j 2 e ctj2 |l - exp(-j 4 tm x )| 



n— 1 n— 1 

< sup sup m' sup sup m m 2 . 



n>l trfT 



n>l ^[H 



i=i 

Hence for d = 1, sup n>1 sup tg r 0T i sup a . g g |t> n (t,x) — t> n,m (t,x)| < Cm"^, and an easy ar- 
gument shows that this inequality can be extended to any d > 1. Furthermore, for any 
m > 1 and t G [0,T], sup n >! swp xeQ E(|w n (t, x) - w n ' m (t, x)\ 2p ) < C ELi^(0> where #i(*) 
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and B±(t) are similar to Bi(t) and B±(t) in the proof of (13. 14p replacing ip(s,y,u(s,y)) — 
ip_(s,K n (y),u(s,K, n (y))) by ip(s, K n (y), u n (s, K n (y))) - ip(A m (s), n n (y), u n (s, K n (y))), B 2 (t) and 
B 5 (t) are similar to B 2 (t) and B 5 (t) replacing <p(s, K n (y),u(s, n n (y))) - ip(s, K n (y),u n (s, K n (y))) 
by <^(A m (s), K n (y), w n (A m (s), - ip(A m (s), n n (y),u n ' m (A m (s), K n (y))) with y? = o or 6 re- 

spectively, and finally B 3 (t) and B 6 (t) are similar to -Bs(t) and B e (t) replacing Gd — (Gd) n by 
{Gd) n — {Gd) n,m - The argument is then similar to that used in the proof of Theorem 13.31 The 
inequalities (12.71) . (1A.13I) . (13.11) and (13.41) provide an upper estimate of B\, (14.371) and (13.11) give 
an upper estimate of B 3 so that B 1 (t) + B 3 (t) < C m'^ l ~^ p . On the other hand, flA.14p and 
fl2~7D show that for some A e]0, 1[, 

B 2 (t) < I (t - s)~ x sup sup E(\u n (A m (s), Kn(y)) - u n > m (A m (s), K n (y))\ 2p ) ds . 

JO n>l y£Q 

A similar argument based on fTA~T3l . fl43U|) . fl3~TD (l3~4"j) proves that B 4 (t) + B 6 (t) < Cmr* for 
any \x e]0, 1[ and shows that for some A e]0, 1[, 

B 5 (t) < fit - s)~ x sup sup E(\u n (A m (s), K n {y)) - u n > m (A m (s), ^ n (y))\ 2p ) ds . 

JO n>l y£Q 

Thus, Gronwall's lemma concludes the proof of (13.231) . The rest of the proof of the theorem, 
which is similar to that of Theorem 13.31 is omitted. □ 



4 Refined estimates of differences of Green functions 

This section is devoted to prove some crucial evaluations for the norms of the difference between 
Gd and its space discretizations [Gd) n , {Gd) n ' m or (G^)™; indeed, as shown in the previous sec- 
tion, they provide the speed of convergence of the scheme. We suppose again that these kernels 
are defined in terms of the homogeneous Dirichlet boundary conditions. Simple modifications 
of the proof yield similar estimates for the homogeneous Neumann ones. 

The main ingredient in the proofs will be the so-called Abel's summation method, which 
is a discrete "integration-by-parts" formula and is classically used in analysis to evaluate non 
absolutely convergent series. More precisely : 

Let (a n ) ne N, (fen)neN be sequences of real numbers, A_i = and A n = J2k=o a k 
n > 0. Then, for any < No < N , one has 

N N N-1 

^ a k b k = (A k - A k _i)b k = A N b N - A No ^b No - ^ A k (b k - b k+1 ) . 

k =N k=N k=N 



In particular, this technique will be employed repeatedly throughout the proofs with x G]0, 2[ 
and a k = cos(/c7ra), for which the corresponding sequence A k satisfies the property \A k \ < 
I . f~sTT , or a k = sin(kTTx), for which A k satisfies a similar inequality, and various monotonous 

| Sm (~2") I 

sequences (bk); see [16J pages 17 - 18 for a more detailed account on the subject. 
Lemma 4.1 There exists positive constants c, C, \x such that for t > 0, n > 2: 

sup \\G d (t, x, .) - (G d ) n (t, x, 0||i < On' 1 (1 + e~ ct , (4.1) 

x&Q 
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r+oo n 

/ sup / \G d (t,x,y)-(G d ) n (t,x,y)\dydt<Cn- 1 
Jo xeQ Jo 



+00 

sup II G d {t,x, •) - {G d ) n {t,x, .)||? a) dt < CrT^ . 
xeQ 



(4.2) 
(4.3) 



Proof : Let 7 > to be fixed later on; the inequalities f]A.9[) . (lA.ip and (1A.11I) imply that for 
< A < 1, 



sup\\G d (t,x,.) - (G d ) n (t,x,.)\\ l( lt < Cn- 2+x , 

xGQ 
p-/n~ 2 

/ sup || G d (t, x, .) - {G d ) n (t, x, .) ||L dt < C n- 2+a 

Jo x&Q 



(4.4) 
(4.5) 



To estimate J" n °^ 2 SU P X 11^(^5 ^5 •) — (G d ) n (t, x, .)|| eft, where || || denotes either the || ||i or || ||( a ) 
norm, we first deal with the case d = 1 and a < 1. 

Case d = 1 and a < 1. As in Gyongy [7J, write |G(t, x,?/) — G n (t, x, y)\ < J2t=iTi(t,x,y), 
where 



Ti(t,x,y) = 








T 2 (t,x,y) = 


E 




Kj<ti-1 


T 3 (t,x,y) = 


E 




Kj<n-1 


T 4 {t,x,y) = 


E 



<^-(x) <pj(y) 



(4.6) 

l<i<n-l 

To study the || ||( a ) norm of a non-negative function R(t, x, .), for x G [0, 1] and i G {1, 2, 3}, let 
A^(x) = {y G [0, 1] : \y — x\ < i n~ or y + x < i rT l or 2 — x — y < i rT 1 } . (4.7) 
Then dy(A l n (x)) < CrT 1 and for x G [0, 1], y,z E A l n (x), \y — z\ < 2in~~ 1 ; furthermore, 
\\R(t, x, .) || J o) < 2 [ || R(t, x, .) 1^)0) 11(a) + \\R{t, x, .) U» (,,-(.) \\\ a) ] . 



Set A^\x) = {{y,z) G Q 2 : |y-x|v|z-x| < 2n- 1 }, ^(x) = {(j/,z) G Q 2 : |y-x|V(x + z) < 



11 /l(2), 



2T2- 1 }, and X 3) (x) = {(y,z) e Q 2 : \y - x\V (2 - x - z) < 2rT x } and for % = 1,2,3, let 



R^\t,x) 



R(t, x, y)\y — z\ a R(t,x, z) dydz . 



(4.8) 



Then II^^Ol^wOII^^CEL^^^)- Let B)t\x) = {{y, z) G Q 2 : 2n- 1 <|j/-x|A 
|z — x\ , 1 2/ — 2! I < 2n -1 } and (x) = {(y, z) G Q 2 : 2n _1 < |j/ — z| A \y — x| A \z — x\} and 
for i — 1 , 2 set 



x 



2 



i?(£, x, — z| a R(t,x, z) dydz . 



(4.9) 



Then \\R(t, x, .)1^2( x ) C (.)|| 2 a j < CY^ i=1 R®(t, x). These notations will be used repeatedly 
throughout the proof for various functions R. 



14 



Estimate of T 2 This term is the most delicate to handle. Set A"(t) := e A ?* — e J 2 "" 2 *. Then 
for any A G [0,2] we have 

< A](t) <C(j/n) 2 fte- cj2t < C n~ A j A e~ cj2 \ (4.10) 
so that flA~7l) with K — A yields 

sup T 2 (t,x,y) < Cn- A r^e- Ct . (4.11) 
xe[o,i] 

Furthermore, T 2 (t,x,y) < \ Y^=i eXj * ¥>j i x ) (y) I + 1 Sj=i e~ j27r2 ' ^ (a;) ^ (y) | , and Abel's sum- 
mation method yields that for 1 < Ni(n) < N 2 (n) < n — 1, 

N 2 (n) 

|.Z^ AJ^J^Kyj L|sin(^)| |sin(^)|J 

Hence for A g]0, 2] and A g]0, -£-r A 1[, we have for any t > 



A+l 
1 



sup / T 2 (t,x,y)dy<Ce- ct n- Ax r A i lx . (4.13) 
xe[o,i] -/o 

In order to bound the || ||( a ) norm of T 2 (t,x, .) for t>^n " 2 , let JVi(n) = N 2 (n) = [n/2] 

and N 3 (n) =n-l. Then T 2 (t,x,y) < Y% =1 T 2 ,i(t, x, y) where T 2tl (t,x,y) = J2 [ j=i and 
for i = 2,3, T 24 (t,x,y) = | EjSl^+Aj (*Vi(*Vi(2/)|- The inequalities (S3Dj) with A = 2 

and (HU) with /? = yield sup^ e[0il] T 2;1 {t,x,y) < C Y}^S n- l e~ cj2t < C n' 1 [1 + H] e" ci . 
Furthermore, sup{T 2 i(t, x, y) ; (t,x,y) G]0, +oo[x [0, l] 2 } < Cy/n. Hence both estimates yield 

sup \\T 2il (t,x,.)\\ 2 (a) < Ce- ct (l + r 1+ f)n- 2+a . (4.14) 

aie[0,l] 

For A e]a, 1[ and /i g]0, 1 - A[, using (14.121) and (14.111) with A = 0, we have for t > 7?i~ 2 , 
with the notations defined in (@~8]) and (S3]): swp xe[0>1] T^(t,x) < Cn a e- ctn2 [l + (nt3)-( A +^)]. 
Similar computations for integrals over the sets An\x) for i = 2, 3 yield sup^Q^ ||T 2j3 (t,x, .) 
xl A 2 iix) (.)\\ 2 a) < Cn a e~ ctn2 . Furthermore, fQ2]) implies that T^t, x) < Ce~ cn ' 2t ( ^ u~ 2 du) 

x ( Jo" 1 v ~ adv ) ^ Ce" cn2 *n a . For (y, z) G £i 2) (x), let I{y, z) < M(y,x) < S{y,z) denote the 

ordered values of \x-y |, and \x-z\. Then f 2 (2 3 } (t, x) < Ce ~ cn2t ! 2n -i< I{ y, z )<s{y,z)<2 [Hv^) 

xM(y, z) j <i?/<iz < Ce _cn The previous inequalities on T 2i 3 yield that for t > 'jn^ 2 , 

sup ||T 2i3 (t,x,.)|| 2 Q) <Cn a e- ct ™ 2 . (4.15) 

x6[0,l] 

We now estimate T 2)2 . Let Co > be a "large" constant to be chosen later on, and suppose 
that t > C n- 2 . For fixed n, t and j G [[y^ + 1, [n/2]], set 0(j) := A"(t). Then 0'(i) = 
2j7r 2 texp[-j 2 7r 2 t] (1 - V©), where for ^ < w < f one sets ^(u) := exp [4n 2 t (w 2 

— sin 2 u) ] . Hence, to apply Abel's summation method, we have to compare ip(u) and 1. 
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Using Taylor's expansion of the functions sine and exponential, we deduce that there exists a 
positive constant C\ such that if C\ = (^ 1 ) 2 , for Co large enough, the map j i— > 0(j) decreases 

on [[y/n\, [n/2]] for t > Let t G [^, ^]. Then there exists a constant C 2 G]0,Ci[ such 

that j h-> 0(j) increases on [[y/n\, [|%]] and decreases on [[|%] + 1, [§ ]] . For t G 

T 2>2 (t,x,y) < Y. 2 i= i T 2,2,i(t,x,y), where one set Bf = ^ anc l 

r 2 ,2,i(*,^y)=| E A i(*)^(^)^-(2/)| + | E KsWnWniv) 

[B-UVi] 

T 2)2 , 2 (t,x, y )=| Xf( t )<Pi(x)<Pi(y) 

j=[B 2 /Vt] 

There exists a constant C > such that if 7 > Co, for every £ > \ 



sup T 2 .2,i(t,x,y) <Cne 
(x,j/)e[o,i] 2 



(4.16) 



For i > let T 2i2 ,i(t, x, y) = T 2>2 (t, a;, y). Using fT4~T2D we deduce that for i > and /3 G [0, 1] 

1 1 



T 2 ,2,i(t,x,y) <C 

<C 

rc Cii 



sin 



• 7v(x—y) y 



+ 



sin 



+ 



7r(x+y) - 
2 • 

1 



.•1 vf J 



i=l 



|x — y| x + y 2 — x — y 



-ctn 



n 



, -2/^-/3-. 

6 ± C n- 2 <i<Cin-l} 



(4.17) 



For i G ^], it remains to bound directly the sum T 2j2)2 (t, x, y). The inequality f!4.1Up 
implies that for B 2 t~^ <j< Bit"*, A](t) < Cn~ 2A t~ A e~ ctj2 for any A G [0, 1]. Therefore, 
the inequality ( 1A.7I) implies that for A G [0, 1], 



sup T 2)2 , 2 (t,x,y) < Cn- 2A r^K 
a:e[o,i] 



(4.18) 



Finally, if Co is large enough, for t G [%, — ] the function ip is increasing on the interval 
[»w/t ' nv^] anc ^ su p{0( w ) : M G [^^, ^j]} < CrT 2 t~ l . Hence Abel's summation method 
implies that for t G [^, ^], 



2 



sup T 2>2:2 (t,x,y) <Cn t 
see [0,1] 



sin 7T 



x—y\ 



+ 



sm 7T 



(4.19) 



The inequalities (I4.17P applied with /3 = | and /3 = 1 respectively and (14.161) imply that for 

, v = A+/i, there exists a constant C > such that for every t > ■yn~ 2 : 



A G]0,a[and/i G]0, 1-A 
su Pxe[o,i] 



T 2 %(t,x)<C 



n -l+a e ~ctn + 1 



{Con-^t^Cin- 1 } 



n 



-4+a+2i/^— ^ 



. Similar computations for 



the integrals over the sets A n (x), i = 2, 3 imply that the same upper estimates hold for 



T£l tl (t,x). Furthermore, (HTTP with (3 g]|,1[ yields Ej=i su Pxe[o,i] T 2 2,i(^ ^) < C( 



n 



-2+a 



xe 



-ctn 



+ n _ 4/3+ a t -2/3 1 



]a, 1[ and /3 gU, 1[, 



{Con-^t^Cin- 1 } 



11) . The inequalities on T^\^ and T^i yield that for z/ G 



0) 



sup ||T 2i2il (t,x, .)||f a) < C[n- 1+a e~ ctn + ! {%<<< a } (r"n- 4+a+2 *' 



ie[o,i 
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For t G ^] OHH]) and (147121 yield that for Ae [0,1], A e]0, a[, [i G]0, 1[ and 1/ = A + //, 
supj-grg^i Y^h=i ^2 ,2,2 x ) — CVi~( 2A+1 ) I ' +a i~( A+ ^ 1 '. Proceeding as for the estimates of T^] we 
deduce that for A = 1 and A g]§, §[ (resp. for < A < §), sup xe[01] T^t, x) < Cn~ 6X+a t- 3X 

and sup s6[0jl] T 2 %(t,x) < Cn~ 4 r 2 - A ( J^-i M ~ 1+X ~ f ^ < Cn-^r 2 . 

The upper estimates of T 1 ^ 2 (£, x), 2 (£, x) and (14.201) imply that for 7 large enough, \i G]|, |[ 
there exists C > such that for every t > ^rT 2 : 

sup ||T 2 . 2 || 2 , < Cn Q \e~ ctn n- 1 + l f c a<t< c lx {n^r* + r 2 n~ 4 + r A n- 2A )] . (4.21) 



xe[o,i] 

Thus (14.141) . (14.151) and (14.211) yield for A G]l, §[, v G]a, 1[ and t > ^n~ 2 with 7 large enough, 



sup \\T 2 (t,x,.W a) <Cn a e" ctn2 +n- 2 e- ct (l + r 1+ f) +n- 1 e- ci " 
xe[o,i] L 

+ 1 { %< i <£i } (n-^^ + n- 4 t- 2 + n- 2X r x 



(4.22) 



Estimates of T 3 Using ( 1A.7I) we deduce that T 3 (t, x,y) <Cn 1 [l + £ x ] e c * . Furthermore, 

set A(l) := [0, 1] n ([^, ^] U [0, ^] U [(2 - A 1, l]). Then dx(A(/)) ^Cr 1 . The 
study of the monotonicity of the function H defined by H(z) = z exp [— 4n 2 tsin 2 (§^)] and 
Abel's summation method yield for large enough 7, t > 7n™ 2 , < A < 1 and y G A(Z), 
sup^i i+i] T 3 (t, x, j/) < C (1 + t 2-)n~ x e~ ct , while for y £ A(Z), sup xg[ i T 3 (t, x, y) < 

f (1 + t~5) e ~ c * [ |y - ^-p 1 + \y + ^-p 1 + |2n - y- ^-p 1 ] • Then, using the partition 
{A(l),A(l) c } and the three previous inequalities we deduce that for A g]0, 1[ and t > r yn~ 2 one 
has 

sup T 3 (t, x, y) < Ce~ ct (l + r^rT 1 . (4.23) 

xe[o,i] 

Similarly, for t > / yn~ 2 and A g]0, 1[ there exists a constant C > such that for every / G 
{0,--- ,n-l}andxe [£, ^], ||T 3 (t, x, .)1 A(0 (.) || 2 a) < Ce~ rf (l + r 1 - A )^ 2+a " 2A . Furthermore, 
when i > 7n~ 2 separate estimates in the cases y,z G" A(Z) and either |y — z\ < n^ 1 or 
|y — z\ > yield that given v g]0,^[, there exists a constant C > such that for every 

Z G {0,... ,n- 1} and x G one has ||T 3 (t, x, .)U(/)4) ||( a ) < Cn^ +Q e - rf (l + r 2 ^). 

These inequalities imply that for t > r yn~ 2 and A g]0, ^[, 

sup ||T 3 (t,x,.)|| 2 a) <Ce~ ct n- 2+ ^ 2X {l+r l - x ). (4.24) 

S6[0,l] 

Estimates of T 4 We suppose that x = -, 1 < Z < n — 1. The general case is easily deduced by 
linear interpolation. For - < y < < k < n — 1, one has k n (y) = - and using ( 1A.7[) . we 
deduce T A (t,x,y) < Cn^^ + l ] e~ ct . Let B{1) := {u G [0,1]; |£ - u\ <\ or ^-u < \ or 2- 
^ - u < as usual, dx{B{l)) < en' 1 . Let then C 1 ^) := {y G [0, l];3u G PI [A; n (y),y]} 
and for i = 1,2, let C*(Z) := {z G [0,1]; 3y G C 1 ^), |y - z\ < i}. Then tfe(C*(Z)) < Cr 1 
and for y ^ C ,1 (0> one nas |y — x| A (y + x) A (2 — x — y) > rT 1 . Computations similar to 
that made to estimate T 3 yield for A g]0, 1[ the existence of a constant C > such that for 
every I G {0, • • • ,n}, and y G C7 2 (Z), Tt(t,l/n,y) < Cn- X (l + r^) e - rf and for y G C 2 (l) c , 
T 4 (t, Z/n, y) < Cn _1 (l + t _ 5) e ~ c * [ |y — ^| : +(y + ^) 1 + |2 — y — ^| X ] . An argument similar 
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to that proving ( 14. 23ft and ( 14.241) implies that for A e]0, 1[ and v e]0, f [, there exists a constant 
C > such that for any £ > 7 n -2 



sup T±(t,x,y) <Ce- ct (l + r — )n-\ 
xe[o,i] 



sup ||T 4 (t,a;, N| " 2 



x6[0,l] 



(«) 



< n -2+a e -ct^ 1 + t -(l+0) n -2i/ - 



Estimate of Ti(t, 2, .). Using (1 A. 71) with /3 = and Jo = n i we have 

sup T x (t,x,y) < C J2 e ~ Ctj2 < Ce~ ctn2 [1 + H]. 
x,»e[o,i] 



(4.25) 
(4.26) 

(4.27) 



On the other hand, since j — > e ■ j2 ' r2 * decreases, Abel's summation method yields Ti(t,x,y) < 
C e~ cn * f . . . 1 X - V „ + . . ; ^c+^m ! • Thus, for AG [0, If and it > we have 

Usmir-/ sin(7r— ■P'-lN L L 



sup / T 1 (t,x,y)dy < Ce^ 2 '[l + f 

see [0,1] Jo 



(4.28) 



Given A e]a, If, \i e]0, If and t > jn 2 , computations similar to the previous ones using the 
partition {A 2 (x), A 2 (x) c } yield that for t > 7^~ 2 , we have 



sup \\T l {t 1 x 1 .)f (a) <Cn a e- cn2t 



2E[0,1] 



(4.29) 



The inequalities (jQgjl . (l4~TBD with A = 2, A = f, and (T 3 + T^)(t,x,y) < CrT l (l + r x )e 
imply the existence of c, C > 0, A e]0, If such that for £ > 0, 



l\„-ct 



sup ||G(t,z,.) -G n (t,x,.)||i < Cr^ 1 (l + r 1 ) e~ ct + r * e" 
xe[o,i] L 



(4.30) 



which implies fO]) for d = 1 with fx = 1. The inequalities (ET38"]) . (I4TT31 with A = 2 and A = |, 
( 14.231) and ( 14.251) with A = \ imply that for some \i e]0, \ [ there exists a constant C > such 
that for every t > 7 n -2 with 7 > large enough, one has 



sup \\G(t,x : .)-G n (t,x,.)\\ 1 <C 

ase [0,1] 



( 1 + ) +e 



-ctn 



l + r 



(4.31) 



On the other hand, the inequalities (T4T291) . (14221) . (I424j) and (OBjl imply that for 1/ e]0, f [, 
A e]l, |[, \x e]a, If, there exist positive constants C and Ci such that for t > jn~ 2 with 7 > 
large enough, one has 



sup ||G(t,x,.)-G"(i,a;,.)|| 2 Q) <CV 

O56[0,l] 



-ctn -ct , f-(\+v) 



n n z 



7? 



2// 



+ ( n- 4+2 ^r^ + n~H~ 2 + n- 2A r A ) 1 



{7«~ 2 <<<Cin- 1 } 



(4.32) 



which proves (14 .11) - (14 .3D for d — 1. 

The general case We extend the inequalities (14.301) . (I4.3ip and (14.321) to any dimension d. 
We use the fact that for any d > 2, we have \Gd(t, x, y) — (Gd) n (t, x,y)\ < J2i=i n», where 

i-1 d 

Ii l =[\[\G{t,x v y j )\)\G{t,x t ,y t )-G n {t,x u y l )\^ J] |G"(t,^-, % )|). (4.33) 

j=l j=i+l 
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Hence, the inequalities (J238), (TATqT) . fj4~30l) . fjOTjl and (14T33]) imply flU) with some // > 0, and 
that for any A e]0, 1[ and any v e]0, l/4[, there exists C > such that for i > 7n~ 2 , 

sup ||G d (t, x, .) - (G d ) n (t, x, .) ||i < Cr v [(l + Hjn" 1 + e~ ctn \\ + *"")] . (4.34) 
x-eQ 

Using (14.41) and integrating (I4.34p with respect to t on [ 7 n _2 ,+oo[ we obtain (14.21) . Finally, 
for 1 < k < d - 1 set a k = a2~ k and set a d = then using (14331 . (P , flA~Tj) . (lA~TTj) and 
(T02D, we deduce that for a e]0,2[, Ci > 0, A e]a,l[, // e]l, §[, v e]0, 24 [ and for i > 7 n" 2 for 
7 > large enough 

d 

SUp \\G d (t,X, . - (G d ) n (t,X, .)|| (a) < C VH 
x-eQ 



x - G n (i, Xj, .)\\( ai ) n*-i 



i=i 



< Cn° 



-ctn z 



+ n -i e - rf " + n" 2 e- rf (r ^ n 2v + r < 1+A ) n~ 2X + l) 
H (^^ ^ + n- A+2 » r 2 +» + n-^ r") l {7n - 2 <,<c in - 1} 1 • (4.35) 



Integrating on [Jj, +oo[ and using (I4.5p . we deduce (14.31) . □ 
We now estimate the norm of the difference {G d ) n and (G d ) n,m . 

Lemma 4.2 Given any T > and v > i/iere exists C > snc/i t/iat 
sup / T / [KW.^^-lG.r^ + Tm- 1 ,!,!/)! 

x&Q Jo J Q 

+ | (G d ) n (t, x, y) - (G d ) n m (t + Tm- 1 ,x,y)\] dydt < C m" 1 ^ , (4.36) 

sup / [\\(G d r(t,x,-)-(G d r> m (t + Tm- 1 ,x,-)\\ 2 {a) 
xgq Jo 

+ \\(G d T(t,x,.) - (G d r m (t,x,.)\\ 2 (a) ]dt < Crrr 1+ * . (4.37) 

Proof : We only prove these inequalities for G n — G n ' m and we at first suppose that d = 1. 
Let Q n ' m = G n - G n ' m and G n ' m = G n > m - G n ' m where G n ' m is defined by (jXlgjl . and C? n ' m is 
defined by 

(nAVm)-l 

^ m (t,x,y)= ]T c^V^ViWy))- ( 4 -38) 

3=1 

Then (IA.20h and (IA.21I) provide upper estimates of the norms of Q n,m . Similar computations 
prove that the same upper estimates hold for the norms of Q n ' m y i.e., for A e]0, |[ and (3 e]a, 1[, 

sup ||^ n ' m (t,x,.)||i <C{l + r x )e- ctm , sup ||^' m (t,x,.)|| 2 Q) <Cm s J L t-te- ctm . (4.39) 
ig[o,i] xe[o,i] 

Let c be a positive constant to be fixed later on; for t < cTm" 1 we estimate separately 
the norms of G n,m (t, x, .) and Q n,m (t, x, .). The inequalities (IA.16I) and (|A.17[) provide the 
estimates of G n ' m . For Q n > m ) we proceed in a similar way. Indeed, j — *> exp(A"t) is decreasing, 
exp(A"t) < e _c * 3 ' 2 for c > 0, and |£ n,m (t, x, y)| < C (n A y'm). Hence the arguments used in 
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the proof of Lemma IA.5I yield that for any c > there exists a constant C > such that for 



t elo. 



cTl 



sup ||^ n ' m (t,x,.) -G n ' m (t,a;,.)||i < C (1 + t" A ) , 

x6[0,l] 

sup ||^ m (t,x,.)-G"' m (t,x,.)||( Q) <C{nAy/m) a . 

xe[o,i] 



(4.40) 
(4.41) 



Furthermore, if * G [cTm' 1 ,! 1 ], |£ n ' m (t, x, y) - G n ' m (t, x, y)\ < f(t,x,y) = f x (t,x,y) + 
f 2 (t,x,y), where ^(t.x.y) = | £ A n,m(t)^(x)^My))\, A^ m (t) = exp m+V>t T 

3=1 ^ 

- (1 - X] and A 2 im (t,x) = exp (\]t) - exp ^ ilxLttel ^ . Using Abel's sum- 

mation method, we have for z = 1, 2, t > — for c large enough, x = - and K n (y) = -, 



\fi(t,x,y)\ < C 



nit 



1 + ^T + 1 



k-t n (j/)| 2—x—Kn{y) 
m 



Furthermore, |T(t,x,y)| < C{n A \/m). 



These inequalities yield that for t > ^ for large enough c and A g]0, 1[, 

||f(t,x,.)||i < cr 1+x m- 1+ ^. 



(4.42) 



For fi G]a, 1[, v g]0, 1 — fi[ and /3 = /i + v E]a, 1[, using the sets /^(^) for i < 3 and 
$ * ^(x) for 7 = 1, 2 and the fact that - < rn= , we deduce that given c large enough, there 

raAvm v ' n — nA^/m' & ° ° ' 

exists constants c, C > such that for every t G [— , T\. 



\\f(t, x, .) || 2 a) < C (n A ^fm) a m~ 2 t~ 2 < C m" 2 +f t~ 2 . 



(4.43) 



For d = 1, the inequalities (TA~20l) . (14T401 and (1Q2l) imply the existence of A g]0, ~[ and 

positive constants c, C such that for any t G]0, T]: 



sup || {{G d ) n - (G d p m ) (t, x, .) ||! < C (1 + r A ) e" ctm + t-i+A^-i 



m 2 



(4.44) 



while the inequalities ijDI) . fl4~39l) . (OB and d4~43l yield the existence of /3 g]0, a A tZ[ and 
positive constants c, c and C such that for every t G]0, T]: 

sup || ((G d ) n - (G d ) n ' m ) (t, x, .)|| 2 < Cm* \{1 + (tm)-^)e^ m + l [£Tm -^ T] (t)m- 2 r 2 } . (4.45) 

Let otk = a 2 k for 1 < k < d — 1 and = at^-i- The inequalities (I4.44|) for d = 1, (IA.9I) and 
(|A~T3|) yield (jOl) for any d, while ( 1051) for d = 1, (lA~Tll and ([Oil yield P5| for any d. 
Integrating with respect to t we deduce the inequalities (I4.36P and (14.371) . □ 



5 Some numerical results 

In order to study the influence of the correlation coefficient a of the Gaussian noise on the 
speed of convergence, we have implemented in C the implicit discretization scheme u n ' m in the 
case of homogeneous boundary conditions in dimension d = 1 for the equation (12.161) . 

To check the influence of the time mesh, we have fixed the space mesh n~ l with n = 500 
and taken the smallest time mesh ttiq 1 with mo = 20736. Using one trajectory of the noise F, 
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we have approximated by the Monte-Carlo method e(mj) = E(|w n ' mo (l, .5) — u n ' mi (l, -5)| 2 ) and 
e(mj) = supj-gro^i E(|n n ' m °(l, x) — u n ' mi (l, x)\ 2 ) for 13 divisors m; of m , ranging from mi = 854 
to ?Tii3 — 144. These simulations have been done for various values of cc, including the case of 
the space-time white noise. Assuming that u n ' m ° is close to u, according to (13.101) and (13.231) . 

these errors should behave like C [m i 2 + n~( 2 ~ a ^] ~ m i 2 for this choice of n and m.j. 
Thus, we have computed the linear regression coefficients c(t) and d(t) (resp. c(i) and d(t)) of 
ln(e(mj)) (resp. of ln(e(mj))), i.e., of the approximation of ln(e(mj)) by c{t) ln(mj) + d{t) as 
well as the corresponding standard deviation sd (resp. sd) for K = 3200 Monte-Carlo iterations 
in the case cr(x) = 0.2 x + 1 and b(x) = x + 2. 

The study of the influence of the space mesh is done in a similar way; we fix the time mesh 
m~ l with m = 32000 and let the smallest space mesh no = 432. Again for various divisors of no, 
using one trajectory of the noise F we have approximated e(rii) = K(\u n ° ,rn (l, .5)— u ni ' m (l, -5)| 2 ) 
and e(rii) = E(|u n °' m (l, .5) — u ni ' m (l, .5)| 2 ) for the 7 divisors n^ of n ranging from 72 to 12. 
Assuming that u n °' m is close to u, according to (13. lOj) and (I3.23p . these errors should behave like 

C [m _ ^ 1_ f ) +n^ 2_Q ' ) ] ~ nj^^ for this choice of rii and m. Thus, we have computed the linear 
regression coefficients -y(x) and 6(x) (resp. j(t) and 6(t)) of hx(e(rii)) (resp. of ln(e(nj))), i.e., 
of the approximation of ln(e(nj)) by j(x) ln(nj) + S(x) as well as the corresponding standard 
deviation SD (resp. SD) for K = 3200 iterations in the case a(x) = 1 and b(x) = 2a; + 3. Both 
sets of results are summarized as follows. 



a 


Theoretical 


c(t) i 


sd 


m 


sd 




exponent 


X = 2 


X = 2 


sup 

X 


sup 

X 


White noise 


0.5 


0.6665 


0.0063 


0.6330 


0.0108 


0.9 


0.55 


0.6954 


0.0121 


0.6853 


0.0130 


0.8 


0.6 


0.7548 


0.0098 


0.7203 


0.0134 


0.7 


0.65 


0.7512 


0.0089 


0.7508 


0.0186 


0.6 


0.7 


0.8158 


0.0143 


0.8007 


0.0090 


0.5 


0.75 


0.8826 


0.0144 


0.8512 


0.0089 


0.4 


0.8 


0.8987 


0.0100 


0.9112 


0.0113 


0.3 


0.85 


0.9592 


0.0117 


0.9135 


0.0117 


0.2 


0.9 


0.9891 


0.0116 


0.9563 


0.0147 


0.1 


0.95 


1.1797 


0.0114 


1.0219 


0.0120 



a 


Theoretical 


j(x) 


SD 


j(x) 


SD 




exponent 


X = 2 


X = 2 


sup 

X 


sup 

X 


White noise 


1.0 


1.2513 


0.0346 


1.2504 


0.0268 


0.9 


1.1 


1.3467 


0.0340 


1.3361 


0.0201 


0.8 


1.2 


1.4347 


0.0336 


1.4251 


0.0211 


0.7 


1.3 


1.5460 


0.0305 


1.5050 


0.0298 


0.6 


1.4 


1.5869 


0.0210 


1.5859 


0.0274 


0.5 


1.5 


1.6714 


0.0280 


1.6671 


0.0272 


0.4 


1.6 


1.7704 


0.0283 


1.7259 


0.0259 


0.3 


1.7 


1.8381 


0.0280 


1.7911 


0.0232 


0.2 


1.8 


1.8978 


0.0274 


1.8503 


0.0208 


0.1 


1.9 


1.9236 


0.0208 


1.9054 


0.0229 



Finally, since our method applies in the case of non-linear coefficients, we have performed 
similar computations for e(mj), e(nij), e(rij) and e(rij) for 1 < i < 13 and 1 < j < 7 with 
K = 3000 iterations in the case u(x) = b(x) = 1 + 0.2 cos(x). The corresponding results are 
summarized as follows 
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a 


Theoretical 


c(t) 


sd 


m 


sd 




exponent 


X 2 


X 2 


sup 

X 


sup 

X 


White noise 


0.5 


0.4915 


0.0602 


0.5200 


0.0431 


0.8 


0.6 


0.5550 


0.0449 


0.6070 


0.0496 


0.5 


0.75 


0.7244 


0.0176 


0.7947 


0.0431 


0.2 


0.9 


0.8607 


0.0225 


0.8571 


0.0429 


a 


Theoretical 


7(x) 


SD 


j(x) 


SD 




exponent 


X g 


X 2 


sup 

X 


sup 

X 


White noise 


1.0 


1.0278 


0.0790 


0.8263 


0.1056 


0.8 


1.2 


1.3628 


0.0830 


1.1276 


0.0684 


0.5 


1.5 


1.5626 


0.0710 


1.5507 


0.0686 


0.2 


1.8 


1.7351 


0.0708 


1.4875 


0.0768 



In this semi-linear case, the speed of convergence is worse and the precision is less than in 
the previous linear case. 

Acknowledgments: The authors wish to express their gratitude to Olivier Catoni and Jacques 
Portes for their helpful advise when writing the C codes. 



A Appendix. 

We start this section with some results concerning the Green kernel Gd in arbitrary dimension 
d > 1. As in the previous sections, we will suppose that Gd and its discretized versions are 
defined with the homogeneous Dirichlet conditions on 5Q; all the results stated remain true for 
the Neumann ones. 

Lemma A.l Let d > 1 and a G]0,2 A d[. There exists some constant C > depending only 
on a, such that for all x, x' in Q = [0, l] d and < t < t' < T : 

sup||G d (t,y,-)||L <Crf , (A.l) 
yeQ 

+oo 

\G d (t, x, •) - G d (t, x', dt <C \x - xf ~ a , (A.2) 

sup/ \\G d (t' -s,x,-)-G d {t-s,x r )\\l a) ds<C\t' -tl 1 ^ , (A.3) 

x&Q Jo 

ft' 

sup/ ||G d (t , - S ,x,-)||( a) ^<C|t , -t| 1 -t. (A.4) 

x&Q Jt 

Proof: To prove (lA~ll . recall the upper estimate of \Gd\ given in (12. 9p . We remark that 
exp ( — c\x — ypt^ 1 ) < exp (— c \y — z\ 2 1^ 1 ) if \x — y\ > \y — z\, while \y — z\~ a > \x — y\~ a if 



o 



| a; — 2/ 1 — \y — z\. Hence (lA.ljl follows from 

sup\\G d (t } x,-)\\ 2 {a) < Cr d ^ + e~ ci rvr a+d - 1 du^j e'^v^ 1 dvj . 

We now prove (IA.2[) and set x 1 — x + v. Then, for < t < \v\ 2 , we have 

\\G d (t, x, ■) - G d {t, x', •) f (a) < 2 [ || G d (t, x, || G d (t, x', •) fj . 
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The change of variables defined by x — y = \v\ T], x — z — |u| £ and t = \v \ 2 s in the first integral 
(and a similar one with x' instead of x in the second one), combined with ( 1A.1I) . yields 

M 2 „ rl 



\G d (t,x,-)-G d (t,x',-)\\l a) dt<C\v 



2-a 



s- a ds 




x \£-v\~ a e~ c — d£d-q + 



{\t-n\>M} 



^77} < G| 



{\t-v\<H} 

2-a 



V 



On the other hand, if t > \v\ 2 , for every j G {!,■■• ,d}, we use the following well-known es- 



timate (see e.g. gTf): -£-G d {t,x,y) 



_ d+i 



< Ct 2 exp I — c 



and the fact that \G d (t,x,y) — 



G d (t,x',y)\ < £^1(^=1!^,^,^ for 
G := G\. Thus, Taylor's formula and for every i G {1 • • • , d} such that v i 7^ 0, the change of 
variables Xj — yj = Vi rjj, Xj — z 3 - = Vi £j for j < i, x' 3 ; — y 3 = Vi rjj, x'j — Zj = Vi £j for j > i + 1, 
and t — vf s yield 

r+00 ^ f+00 pi 

/ \\G d (t,x r )-G d (t,x' r )\\l ) dt<Cy2l {Vim \v l \ 2 ~ a / s-^ds d\ 



1-7,-r 



if,-r 



' d^jdrjj 



Splitting again the integrals between — rjj\ < \rjj\} and — 77.7 1 > \r)j\} for j ^ i and 
{|& - %| < + A|} and {|& - r/;| > |r/; + A|} yields 



+00 
M 2 



|G d (t J x,.)-G d (t,ar' > .)||J a) <C|v 



2-a 



This completes the proof of ( 1A.2I) . On the other hand, f lA.31) is obtained using similar arguments 
and the change of variables defined by t — s = hr, y — x = y/hr] and z — x = Vh^ (where 



d+2 Jg-1>K 



□ 



h = t' — t > 0), Taylor's formula and the estimate | J^G^t, x,y)\ < Ct~~e 

We recall the following well-known set of estimates. The proofs can be found in [TS] for 
d = 1 and are easily deduced for any d > 2. By convention set Gd(t, x, = if t < 0. 

Lemma A. 2 For x, x' G Q, < t < f < T and // G]0, 1[, 

"+00 /> 

/ \G d (t,x,y) - G d (t,x',y)\dydt < C\x — x'\ 
Jq 

[ [ \G d {t-s,x,y)-G d {t' -s,x,y)\dyds < C \t' - t\^ 
Jo Jo 



(A.5) 
(A.6) 



The following technical results are needed to obtain refined estimates for the discretized 
kernels G n and G n ' m . The proofs, based on simple comparison between series and integrals for 
piecewise monotone functions, are omitted. 

Lemma A. 3 For any c > there exists a constant C > such that, for K > 0, (3 G [0, 1[ ; 
t > 0, a > 1 and Jq > 1, 



1 +t 2 



(A.7) 
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j=Jo 



cTj 



2 \ -i 



ni 



j=Jo 



cTf 



m 



< 



Cm 



1 + 



cTJl 



JoT(a — 1 

The following lemma bounds the || ||i and || || Q norms of (Gd) n (t,x, .). 



m 



(A. 



Lemma A. 4 There exists a constant C > such that for every t > 0, d > 1, X > and 

< a < f3 < d A2, 



sup sup || (G d ) n (t, x, .) || ! <C (1 + r A ) e 

n x£Q 

supsup||(Gf d ) n (*,z,OII(«) <Ct~^e- ct , 

n x£Q 

sup\\{G d ) n (t,x,.)\\ 2 ^ <Cn a e- d 



xeQ 



1(a) 



(A.9) 
(A.10) 
(A.11) 



Proof: It suffices to check these inequalities for x = l/n, < I < n. We at first prove them 



for d = 1. Using Abel's summation method, we have, since j 



X' 



decreasing, for k = K n (y), 

\G n (t,x,y)\ <Ce A "*| sin ( 



-4n 2 sin 2 (f) is 



x ~ n n (y), - 1 

7T 



sin (tt- 



x + K n (y) 



-i 



(A.12) 



Fix < A < \ and t > 0; (TAT9]) and (TAT7jl with /3 = 0, J = 1 yield (TA~9l for d = 1. 
Let be the sets defined by gUD and for < / < n, set = A^l/n). Then 

da; ^L>n ) (/)J < ^ and, if y g D ( n\l), one has \x - k n (y)\ > -\x - y\, \x + k n (y)\ > ~\x + y\ 



and similarly, if y ^ Dn (I), \x — k n (y)\ > -\x — y\, \x + k n (y)\ > ~\x + y\. Thus, for 
every n > 1 and < I < n, \\G n {t, l/n, .)\\ 2 {a) < C(T X + T 2 + T 3 ), where T h 1 < % < 3 
is the integral of \G n {t, l/n, y)\ \y — z\~ a \G(t, l/n, z)\ respectively on the set A\ = {(y,z) : 
y G D®(l),z G D%\l)}, A 2 = {(y,z) : y G D^{lf, z G D^{lf, \y - z\ < n~ 1 } and 
A 3 = {(y,z) : y G D ( n\l) c , z G D { n\l) c , \y - z\ > n' 1 }. Thus, using ( fA~7l) and (fA~T2l . we 
deduce upper estimates of Tj for 1 < z < 3 which imply ( 1A.10I) when d = 1. 
Again, to prove (IA.110 . it suffices to show that sup 0<Kn \\G n {t,l/n, .)\\ 2 a ^ < Ce~ ct n a . Using 
the sets A i: 1 < i < 3, the inequality ( 1A.12I) . the crude estimate \G n (t,x,y)\ < Cne~ ct , and 
replacing in products involving two of the terms \x — y\~ l , \x — z\~ l and \y — z\~ 2 the largest 
norm by the smallest one, a similar computation yields for a < X < 1 and < // < 1, 

\\G n {t,x, .)|| 2 Q) < C [ n A +f n - A -f+« + n a ] e" c * < Cn a e- d . 

Since (G (i r(t,x, Z /) = ntiG rn (i , Xj, yi), ( 1A.9I) for d = 1 immediately yields 1 ) A. 91) for any d. For 
d > 2, and 1 < i < d — 1, set = a2" ! and set = a^-i- Then using (12. 2p . the inequality 
(TA~T0|) (resp. (jATTT]) ) for d = 1, we deduce (IA~T0|) (resp. (jXTTB ) for every d. □ 
We now prove a similar result for the norms of (Gd) n ' m {t,x, .). 

Lemma A. 5 For every X g]0, ^[and (3 G]a,d A 2[, there exist positive constants c and C such 
that for every t g]0, T], 



1 



sup||(G d ) n ' m (t,x,.)||i<Ce- ct (l + r A ) 
sup || (G d ) n ' m (t,x, 0||( a ) <Ce~ ct \(n A \pm) a A (l + t~ 2 ) 



-ctm 



(A.13) 
(A.14) 
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Proof : For m > 1, set (G d ) n > m (t, x, y) := (G d p m (t, x, y) - (G d ) n > m (t, x, y), where 



_ _ _ — _ mt 

(G d r> m (t,x,y) = E II i 1 -Tm^Xl)-^ (A-15) 

ke{l,- ,[(nAVm)-l)]} d i=1 



Let (Gi 



G" 1 '" 1 and (d 



G" 1 '" 1 . Since j -> (l - £ A") [ T 1 is decreasing, |G n ' m (t + 



Tm 1 ,a;, y)\ < C(n A yfm)e c * and Abel's summation method yields that for x 



and 



K n (y) = t- \G m ' n (t + Tm-\x,y)\ < Ce 



ct 



+ 



+ 



. Finally, since for 



\x—Kn(y)\ x+K n (y) 2-x-K n (y) 

j < y/m, In (l + ^4n 2 sin (g) ) > CfTm' 1 . Using itOl we deduce IG^t + Tm- 1 ,!,^ < 

C ^]^ v/ ^~ 1 e" cfj2 < Ce~ ct (l + t~?). Thus, repeating the arguments used to prove ( 1A.9I) - 
( 1A.11I) we deduce that for A > 0, < a < (3 < d A 2, 



sup \\G n ' m {t, x, .) ||i <C e - ct (1 + r A ) 



sup||G"' m (t,x,.)||( a) <Ce- ct 



(l a (n A y/m) a 



(A.16) 
(A.17) 



We finally give an upper estimate of the norms of (G d ) n ' m (t,x, .) and thus we suppose that 
y/m < n. Using (1A.8D we deduce the existence of positive constants c, C such that for 
t < f , sup^ 6Q |G"' m (t,x,y)| < a/^=r( m T)-J(l + y 2 )- ([¥l+1) ^- Hence for t < 2Tm"\ 
since [^] = 1 or 2, for x,y G Q, |G n ' m (t, x, y)| < C < y/m < H while for t > 2Tm~\ 
\G n > m {t,x,y)\ < C(mT)-3 + £/ 2 )~ ([ ^ ]+1) A/ < Cr 1 ^-* < Cri This implies that 



sup \G n ' m {t,x,y)\ < C(l + H) . 

Furthermore, j — > (1 — TA^/m) - ^ decreases and (l — TA^/m) 
foi x = l/n and = k/n by Abel's summation method 

1 1 



(A.18) 



< C e- cTm . Hence 



\G n ' m {t + Tmr l ,x,y)\ < Ce 



ctm 



+ 



\x-K n {y)\ x + K n (y) 



+ 



s - K n (y) 



(A.19) 



An argument similar to that used to prove ( 1A.9I) implies that for A g]0, |[, there exists C > 
such that for t G]0,T[, 

sup||G n ' m (t,x,.)||i < C(l +r A ) . (A.20) 



Finally, for x = l/n let D l m (l) = {z G [0, 1] : |x— z| < iy/m, or x+z < z-y/ro, or 2—x—z < iy/m}. 
Then since n > y/m, for y G" -D^(Z) we deduce that \x — «n(y) > \\x — y|, |x + « n (j/) > ||x — y| 
and \2 — x — K n (y) > ^\x — y\. Hence, the arguments used to prove (lA.lOl) with d = 1 and ( 1A.11I) 



show that \\G n ' m {t+Tm- l ,x, .)1 D3 



(0 c H(a) 



< Ce 



-ctm 



ni 



2 . Furthermore, ( 1A.18I) and ( 1A.19I) imply 



that the same upper estimate holds for \\G n,m (t+T m 1 , x, .)1.S3,(0 IIL)- These inequalities imply 
that for /3 G]a, d A 2[, 



sup ||G n ' m (t,ar,. , )||L<Ce- < "" 



xe[o,i] 



1(a) 



f [l + (*m)~a] . 



(A.21) 



Hence ( 1A.21I) and (1A.17|) imply that (1A.14|) holds for d — 1. Finally, as in the proof of Lemma 
lA~4l flA33|) and flA~14j) hold for any d > 1. □ 

We finally prove upper estimates for the norms of time increments of {G d ) n and set by convention 
(G d ) n (t,x,.) = 0tft<0. 
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Lemma A. 6 For any T > 0, there exists C > such that for any h > 

"t+h 



pt+n 

sup sup sup / \\(G d ) n (t- s,x,.) - (G d ) n (t + h-s,x,.)\\xda < C h* , (A.22) 

n>l x&Q te[0,T] Jo 

rt+h 

supsupsup / \\(G d ) n (t- s,x,.) - (G d ) n {t + h-s,x,.)\\ 2 {a) ds <Ch l ~% . (A.23) 

n>\ x£Q te\0,T] JO 



Proof: Inequality (1A.9I) implies that for x G Q, X > and n > 1, J^ +h J Q \ (Gd) n (t + h — 
s, x, y)\dyds < C L s~ x ds < Ch 1 ~ x . Using ( 1A.11I) we deduce that for any c > and h < 



"£h rch 



n 



pen pen 

/ \\(G d ) n (s,x,.)\\ 2 (a) ds< n a ds<Ch l -%. (A.24) 
</o Jo 

Suppose that s > n~ 2 and that d — 1. Then by (IA. 7[) . |G n (s, x, y)| < C (1 + s~l). Then using 
( 1A.12I) and proceeding as in the proof of ( lA.llj) . replacing the sets A l n {x) defined by (14.71) by 
the sets A\{x) = {y G [0, 1] : \y — x\ < i\/h or y + x < i\fh or 2 — x — y < iVh}, since we 
have assumed that n~ l < y/s, we deduce that ||G n (s, x, .) < C s~2 e~ ct . Let a& = a2~ k 
for 1 < k < d — 1 and a d = a d -i] using the inequality (12. 2L we deduce that for s > n~ 2 , 
then \\(G d ) n (s,x,.)\\ 2 a) < Cs~^e~ ct . Hence this inequality and flA~24l) imply J t t+h \\ (G d ) n (t + 

h — s,x, .) ||? a )rfs < C"[ J Q A " n°<is + J^ An -2 s~?<is] < Ch x ~% , which yields f)A.23[) . To complete 
the proof of ( 1A.22I) and ( 1A.23I) . set t' = t + h and consider the integrals on the interval [0,t[. 

Then for d = 1 and x = l - one has \G n {t,x,y) - G n {t',x,y)\ = £J=i e'^^ te) 
g -4n fism (^) ip^ x }ipj{K n {\)) . Thus (IA.7p implies the existence of C > such that for any 

ri > 1 andx = ln-\ f* J Q \G n {t-s,x,y)-G n {t' -s,x,y)\dyds < C J' 2 [{fh) Al] < Chi, 
which proves (IA.22p . The inequality (1A.24I) proves that given any c > 0, 

sup sup / (\\G n {s ) x ) )\\ 2 {a) + \G n {s + h,x,)\\ 2 {a ^)ds<Ch 1 '^ . (A.25) 

n>l x6[0,l] Jo ^ ' 

Fix c > large enough, let t > ch and set $(j) = exp(— 4n 2 t sin 2 (||)) — exp(— 4n 2 t' sin 2 (||)) > 
0. Then $'(j) = 2nsin (^) [tf exp ( - 4n 2 f sin 2 (g)) - t exp ( - 4n 2 t'sin 2 (£))]• Then the 
arguments used to estimate $2(^) and then x, y)\ in the proof of Lemma [4.21 show that 
there exists C > such that for any s G [ch,T], x = ln~ l and y G [0, 1], | YT~ h -\ {e XjS — 

e^ +h ^(x) <fMn{y))\ < C exp (A»_ 4]S ) + + ^^j] , while « implies 

that IE™! 1 -! (e A " a ~ e A " ( ' +h) )vi(a;)^(K n (y)| < Cs~i exp (A n _j s). Using again the sets 
^(x), we deduce that for any s G [c/t, T], sup n>1 sup^Q^ \\G n (s,x, .) — G n (s + h,x, -)\\ 2 a \ < 

Ch 1 ^^. Thus, sup n>1 sup^gjQ ^ f^ h \\G n (s,x, .) — G n (s + h,x, .)\\ 2 a ^ds < Ch 1 '^ . This inequality 
and ( 1A.25I) yield (1A.23I) for d = 1. We extend the lemma in any dimension d > 1 as in the 
proof of Lemma IA.41 □ 
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